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Abstract Recently, an extension of the geodesic equations of motion using 
the Dark Energy length scale was proposed. Here, we apply this extension 
to the analyzing the motion of test particles at the galactic scale and longer. 

XV . A cosmological check of the extension is made using the observed rotational 

velocity curves and core sizes of 1393 spiral galaxies. We derive the density 

_. . profile of a model galaxy using this extension, and with it, we calculate as to be 

\^ . We then calculate i?2oo to be 206±53 kpc, which is in reasonable agreement 

(-H ' with observations. 

o 

L_i i 1 Introduction 

C^ . In a previous paper [T] , we constructed an extension of the geodesic equations 

of motion (GEOM). This construction is possible because with the discovery of 
Dark Energy, Ade = (7.2lto:84) x IQ-^Og/cm^ 0-0], there is now a length 
^ [ scale, Xde — c/{AdeG)^/'^, associated with the universe. As this length scale 

is also not associated with the mass of any known particle, this extension 
fv^ ', does not violate various statements of the equivalence principle. Importantly, 

V/^ ' the extension does not change the GEOM for massless particles, and thus as- 

tronomical observations of the universe — which are based on the trajectory of 
^^ ' photons — remain unchanged. At 14010120 Mpc, the shear scale of Xde ensures 

f^ I that effects of this extension will not have already been observed either in the 

motion of bodies in the solar system, or in terrestrial experiments. Indeed, by 
analyzing the effects of the extension at these length scales, we established a 
lower bound to a a- This a a is the only free parameter in the theory, and gives 
. , the power law dependence of the extension on the ratio c^R/Ae>eG, where R 

j^ ' is the Ricci scalar. 
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At the conclusion of [T] , we argued that it is only at galactic length scales 
and longer will we expect effects of the extension to be relevant. A study of 
these effects is the focus of this paper. In particular, we analyze the restric- 
tions on our extension of the GEOM due to both observations of galactic 
structure, and recent measurements by WMAP. We find that these observa- 
tions and measurements do not rule out our extension of the GEOM. To the 
contrary, by applying the extension to an analysis of the motion of stars in 
spiral galaxies, we are able to calculate erg — the rms fluctuation in the density 
of matter at 8h~^ Mpc — to be 0.73±o.i2- This value is in excellent agreement 
with the WMAP measurement of 0.761^g q48 f*^^ '^s- ^^ ^^^ ^^^^ ^^^^ ^^ '^^l" 
culate i?200 — the distance from the center of a galaxy at which the density of 
matter equals 200 times that of the critical density — to be 206±53 kpc, which 
is in reasonable agreement with observations. Calculations of both quantities 
are possible because of an unexpected connection between the effects of our 
extension on the motion of test particles at galactic scales with cosmological 
length scales. This connection allows us to set a definite value of 1.56±o.io for 
aA that is equal to the rough lower bound on a^i established in [T] for effects 
of the extension to be unobservable at solar system and terrestrial scales. 

This stringent test of the extended GEOM is only possible because of the 
scale of Xde- With a value of 14010|2o Mpc, it is expected that any effects 
of the extension will only be apparent at the galactic scale or longer. (Indeed, 
Xde is so large that it is only because of the nonlinear dependence of the 
extension on R that effects at the galactic scale are apparent at all.) It is 
also precisely at the galactic length scale where deviations from motion under 
Newtonian gravity appear. Using a simple model of a spiral galaxy, we are able 
to determine the density profile of the galaxy by applying our extension of the 
GEOM to the motion of stars within it. We find that effects of our extension 
necessarily extend beyond the galactic scale. Not unexpectedly, at these scales 
the 1/r interaction potential that is expected between galaxies from Newtonian 
gravity is now logarithmic, which is consistent with the interaction potential 
between galaxies and galactic clusters inferred through observations. At even 
larger length scales, we find that the predicted density profile for the model 
galaxy goes to zero exponentially fast at distances beyond the Hubble length 
scale, Xh — c/H, (where H = hHo, and Hq = 100 km/s/ Mpc) from the 
center of the model galaxy. While this result is certainly physically reasonable, 
it is surprising that such a length scale naturally appears in the theory, even 
though a cosmological model is not mentioned either in its construction, or in 
its analysis. This unexpected connection with cosmology allows us to set the 
value of UA to be 1.56±o.io using the WMAP measurements of h and Ade- 

We have also used this density profile to calculate explicitly as- This cal- 
culation is possible because of four data sets in the literature, [5] - [13] . These 
data sets are the result of observations — made over a thirty-year span — of 
galactic rotation curves that give both the asymptotic rotational velocity and 
the core sizes of 1393 spiral galaxies. While only a very small fraction of the 
observed galaxies in the universe, the size of the data set and the fact that 
the vast majority of these observations were unbiased, allows us to obtain av- 



erage values for the parameters used to characterize the model galaxy. These 
parameters, along with the value of a a, is used to calculate a cts that is in 
excellent agreement with WMAP measurements. We are also able to use these 
parameters to predict a value of 206±53 kpc for i?200i which is in reasonable 
agreement with observations. 

Like Peebles' model for structure formation [Tl], the total density of mat- 
ter for our model galaxy can be written as a sum of a asymptotic, background 
density, Pasymp(x), and a linear perturbation, p}j(x). The matter that makes 
up Pasymp(x) docs uot Contribute to the motion of stars within the galaxy, 
while the matter that makes up p}j{x) does. Unlike Peebles' model, however, 
Pasymp(x) is not a constant, but instead varies inversely with distance from 
the galactic core, and dies off exponentially fast beyond the Hubble scale. 
Moreover, the form of Pasymp(x) depends only on the dimensionality of space- 
time and a A, while the scale at which it decreases depends only on Xde- In 
this sense, /3asymp(x) is universal, and does not depend on the detail structure 
of the galaxy. This is in contrast to p}j(x), which depends explicitly on the 
structure of the galaxy near the galactic core both in form and in scale. 

While /3asymp(x) docs not contribute to the motion of stars within galax- 
ies, it does contribute to the deflection of light. By construction, the extended 
GEOM does not affect the equations of motion for massless particles. Light 
still travels along geodesies, and the degree of the deflection of light is de- 
termined by the total local density of matter. As such, using the deflection 
of light to measure the local matter density will result in a measurement of 
/9asymp(x) -|- yo}j(x). In Contrast, the motion of stars in galaxies are affected 
by our extension of the GEOM. As such, their motion is determined solely 
by p}j(x), and thus when this motion is used to determine the local matter 
density, what is measured is only the matter that makes up p}j{x). As we find 
p\j « Pasymp outsidc of & fcw galactic core radii, the presence of the vast 
majority of matter in the universe can only be inferred by though the effects 
that the local density of matter has on the trajectory of light. 

The rest of the paper is divided into five parts. In the first part, an overview 
of the extended GEOM is given, and the properties of the extension needed 
in this paper is outlined. In the second part, we introduce our model of the 
galaxy, and using the extended GEOM, we derive the density profile of the 
galaxy given a rotational velocity curve for it. We show that this density 
dies off exponentially fast beyond a fixed distance from the galaxy. In the 
third part, we use WMAP measurements of h to determine a a, and calculate 
(78 and i?200 using the density profile of our model galaxy and observational 
data on galactic rotation curves from the literature. The value for erg is then 
compared with WMAP measurements. In the fourth part, we calculate the 
gravitational potential for the model galaxy, and determine which portion of 
the density can be determined through direct observations of the motion of 
stars in the galaxy, and which can only be determined through the deflection 
of light. Concluding remarks can then be found in the last part. 



2 An Overview of the Extended GEOM 

While there is currently no consensus as to the nature of Dark Energy, modi- 
fications to Einstein's equations to include the cosmological constant are both 
well known and minimal. In addition, WMAP measured the ratio of the pres- 
sure to energy density ratio for Dark Energy to be — 0.9671q'q72; this is within 
experimental error of —1, the ratio expected for the cosmological constant. We 
thus identify Dark Energy with the cosmological constant in this paper, and 
require only that Ajje changes so slowly that it can be considered a constant 
in our analysis. Einstein's field equations are then 

R _i p I AdeG _ 8ttG ^ 

where T^i, is the energy-momentum tensor for matter, R^i^ is the Ricci tensor, 
Greek indices run from to 3, and the signature of g^v is (1,-1,— 1,-1). 

The extended GEOM for a test particle with mass, rra, is obtained from 
the Lagrangian 

£Ext = ^yKlc^R/ADEG] {9^,,x^'x'y- . (2) 

While the function IH is arbitrary, in [I] we argued that the simplest choice is 

$n(a;) ^ [l + £>(a;)]^/^ where 



while 



^{-)-x{a^)l YT^TT^' ^^^ 



ds _ sin[7r/(l + ua)] , . 



X{aA) Jo l + si+"-i 7r/{l + aA) ' 

is defined so that D{Q) = 1. Here, a a is a constant, and is the only free 
parameter in the theory. In [T] , we showed that for the effects of the extension 
not to have already been observed in terrestrial experiments, aA must be 
between 1.28 (for Ade = 10"^^ g/cm^) and 1.58 (for Ade = 10"^^ g/cm^). 
While the equations of motion derived from £ext is 

-g^^c'(g'^-~^]wAogy\[c'R/ADEG], (5) 

we are interested in the motion of stars in galaxies. Moreover, WMAP and the 
Supernova Legacy Survey put Qk = — 0.011±o.oi2, and thus the universe is es- 
sentially fiat. Indeed, WMAP's value for h is determined with this assumption. 
As such, we are working in the nonrelativistic and weak gravity limits, and 
therefore take the metric to be g^i, = rj^n -\- h^^. Here, rjp^u is the Minkowski 
metric, and /i^^ is a small perturbation. The only nonzero component of /i^j, 
is /loo = 2^/c^, where <? is the Newtonian gravitational potential. Equation 
(P) then reduces to 

V2<|> + 2^^^f^'? = 47rpG-ylD£;G', (6) 



with the additional terms due to the cosmological constant. 

As we are deaUng with the motion of stars in galaxies, taking 77^^ to be the 
background metric would seem to be straightforward. There is one subtlety, 
however. Even when T/^^^ = in Eq. 1^, Ade is present, and the spacetime 
is not flat; at scales comparable to Xde^ tbe spacetime will be significantly 
different from Minkowski space even in the absence of matter. At I4OIOI20 
Mpc, \de is over three times the Hubble scale, Xh, however, and taking the 
background metric to be flat is a good approximation throughout most of the 
physically relevant length scales. By restricting ourselves to length scales much 
less than Xde , taking rj^^, as the background metric is a good approximation, 
and the terms proportional to Ade in Eq. ([6]) can be neglected. 

Using Eq. (O, Eq. ([5]) reduces to 




l+aA ^ ^ 



Vp, (7) 



in the nonrelativistic and weak gravity limits. Here, we have related R to the 
trace of the energy-momentum tensor, T, using Eq. ([T]). We have also assumed 
that the spacetime is spatially symmetric, and that the particle moves through 
an ambient, nonrelativistic fluid with density, p. We are dealing with only 
gravitational forces, and thus do not differentiate between baryonic matter 
and Dark Matter in p. As shown in [1], the energy-momentum tensor for the 
fluid can be approximated as T^^, sa pu^Uj, in the nonrelativistic and weak 
gravity limits even though elements in the fluid propagate under the extended 
GEOM instead of the GEOM. Because Ade ^ 10~^° g/cm^, we expect that 
Sttp/Ade > in galaxies, and thus have also used the expansion 

D(4 + 8.pM..)^xE „ni:,^.. 4+-£ , (8) 



n=0 



n(l -I- oa) + ayi V ^de) 



in obtaining Eq. ([7]). 



3 Effects on the Galactic Scale 

While definitive, a first principle calculation of the galactic rotation curves 
using the extended GEOM to describe the motion each star in a galaxy is 
analytically intractable. Instead, given a rotation curve for a model galaxy, 
we will use the extended GEOM to derive the density profile for the galaxy, 
and through this profile, describe the motion of test particles within it. Not 
surprisingly, we find that inside the galactic core stars move as though they 
are in a Newtonian potential. Outside of the core, on the other hand, stars 
move as though they are in a non-Newtonian, logarithmic potential. 

Since the lower bound on aA is between 1.28 and 1.58, we take ua = 3/2 
as a guide when making approximations in this section. In the next section, 
we will use WMAP measurements of h and Qa to determine a a- 



3.1 The Model Galaxy 

A number of geometries have been used to model the formation of galaxies |15j 
- [17) . including a spherical geometry. Because we will be making connection 
with cosmology, we are interested in the large-distance properties of the density 
profile, and at such distances, detail structures of galaxies are washed out; only 
the spherically symmetric features survive. We thus use a spherical geometry 
to model our idealized galaxy, and will focus on the general structure of the 
galaxy instead of on the details. 

To determine the density profile, we divide space into the following three 
regions. The first. Region I = {r | r < th, and p ^ AdeI'2-it}i encompasses 
the galactic core with radius, rn- Here, the density of matter is much larger 
than yl£)E/27r. The second, Region ll — {r\rH<r< ru, and p ^ Ade/^tt}, 
encompasses the region outside the core that contains stars undergoing rota- 
tional motion with constant velocity. Here, the density of matter is also greater 
than Ajje/'^'^, and the region extends out to a distance of r//, which is deter- 
mined by the theory. While this region of constant rotational motion is much 
larger than the size of the galaxy, there are no other galaxies in the model, and 
we find that there is a cutoff in the density at distances beyond r//. Finally, 
the third. Region HI = {r \ ru < r, and p ^ Ajje/^t^}, encompasses the 
region where the density of matter is less than Aue/^tt. We will find that the 
density is exponentially small here, and we will see in the next section that 
this region is not physically relevant. 

We assume that all the stars in the model galaxy undergo circular motion. 
While this is an approximation, galactic rotation curves are determined with 
stars that undergo this motion, and we use these curves as inputs for our 
analysis. With this approximation, the acceleration of each star, a = x, is a 
function of r only. We can then take the divergence of Eq. ^, and obtain 






f{r)=p-X>^DE- 7 : \ 1+^^ >, -, (9) 



where /(r) = — V • a/AnG is considered a source. Because 

-V-a=l|:Hrr], (10) 

given a velocity curve for the galaxy, v{r), a can be found, and the source 
determined. 

While a number of different models have been used in the literature to fit 
the observed rotational velocity curves, we use a particularly simple idealiza- 
tion of the curves: 

ideal /„^ _ f vht/th for r <rH 



' vh tor r > Tji 



where vh is the asymptotic velocity of the velocity curve. We make this ide- 
alization because our purpose is not to derive the precise density profile of a 
galaxy; our choice of geometry does not allow for such a derivation. Instead, 
our purpose is to determine the overall features of the density profile; for this 
purpose, -yideaij-y,-! jg gyfgcient. 

While ti"i<=^'(r) is continuous, f{r) is not. We find that 



pij = 'iv'jj / AttGt'Jj for r < th 
pijr^j/ir^ for r > rn- 



/W - ^ ....2 /q^2 



Here, pn is identified as the density of matter in the galactic core, which is a 
constant for our idealized velocity curve. 



3.2 The Density Profile of the Model Galaxy 

In both Regions I and II, p ^ Ade/^tt, and Eq. (jH]) may be approximated as 



where u = r/x^''^^DE, and V„ denotes the derivative with respect to u. 

In Region III, on the other hand, we find that /(r) is negligibly small, and 
can be set to zero. Equation ^ then reduces to 

= p-- TT^V^p. (12) 

' I -\- 41+0,1 "'^ ^ ' 

3.2.1 The Solution in Regions I and II 

Because f{u) = pn in Region I, we find that the density, pi{r), in this region 
is simply pn. 

To find the density, pii{r), in Region II, asymptotic analysis |5D] is used. 
We first find a solution to Eq. (ITT]) in a region where PH(?'///r)^/3 ^ p{t)- 
This is equivalent to making the anzatz that within Region II there is a point 
''asymp bcyoud which the driving term in Eq. (|lip is negligibly small compared 
to the density. Equation dTTT) then reduces to 

where Pasymp is the asymptotic density. We then perturb off Pasymp to find the 
density near the galactic core. 

The solution to Eq. ([13]) has the form 

Pasymp = -3 S{aA)vP , (14) 



where p and S{aA) are solutions to 

p(l - aAp) 



= 1- 



1 



S{aA) 



l+OA yp(l+ayi)+2' 



(15) 



Thus p = —2/(1 + O/i), while S{aA) is given by 

2(1 + 3a^) 1 



0= 1 



(l + a^)2 [l]{aA)] 



l+Q^ 



(16) 



For Pasymp to be positive, S{aA) > 0, and thus there must be positive solutions 
to Eq. (dni). Such solutions exist only if a a is the ratio of a odd integer to an 
even integer. Our choice of ua — 3/2 satisfies this criteria, and we arrive at 
the following asymptotic solution 



Pasyinp V j 



J^DE r-./ N fX^DE 



Sir 



-r(ayi) 






1/il + aA) 



with 



S{aA) = 



2(l + 3aA) 



(1 + «a)^ 



l/{l+aA) 



(17) 



(18) 



To justify our anzatz that rasymp lies within Region II, we set /(^asymp) 
/Oasymp(fasymp), and find the bound 






TH 



ZA 



DE 



S{aA) 



(19) 



where uh — th/x^^'^^de- For a a — 3/2, pn ^ 10^^^ g/cTc?, and rn = 10 
kpc, rasymp < 6.8r//; the anzatz is thus valid throughout the great majority of 
Region II for the range of galaxies we are interested in. The upper limit, r//, 
to Region II, on the other hand, is found by setting 87rpasymp(''//)/^DS = 4. 
This gives 



rii 



l^io^A) 



(l+ayl)/2 



X 



i/^A 



DE- 



(20) 



For a^ = 3/2, m« 0.27 AzjB. 

Next, to see the structural details of the galaxy, we take pn = jOasymp + zo)/- 
Expanding Eq. dill) to first order in p\j gives 



2(1 + 3a^) 
{l + aAY 



Ph / Uh 
3 \ u 



^^3^%+V^.pA. (21) 

(1 + q;a) m 



where p^j — u^pjj- The particular solution to Eq. (PT|) is again the constant 
solution, but now for p/^; this corresponds to p]j ~ puT^jir^ ■ The solution 



Stt \ H 


\ V(i- 
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r"D 


2 
+ 




(^ ^ j 1 CcosCOS 


I'log 


r r \ 


+ Csin sin 


i^log 


f r \ 



of the homogeneous equation is straightforward, and combined with Pasymp 
from Eq. (J17p . we arrive at the density profile in Region II 

Piiir) 

/T'it\^/^/ /t"\ /t"\\ 

(22) 
/J L \''H / j/ 

Here i^ = [2(1 + 3ayi)/(l + ayi)^ - 1/4] ^^^ while 

2 

Ccos = 7:PH — Pasymp (?■_?/), 

7 l(l + 5aA) , . , . 

l^Uin = i;PH - X— — ^Pasymp(r-//), (23) 

6 2(1 + aA) 

are determined by the boundary conditions for p{r) at rn- 

The density, pii{r), thus consists of the sum of two parts. The first part, 
/3asymp(?'), corrcspouds to a background, asymptotic density, and depends 
solely on Dark Energy, fundamental constants, the exponent a a, and the di- 
mensionality of spacetime. It is universal, and has the same form irrespective 
of the detailed structure of the galaxy. The second part, p]j{r), on the other 
hand, does depend on the detail structure of the galaxy. Variations in p]j are 
on the scale of r^, the core size; in contrast, variations in Pasymp are on the 
scale of XjjE, the Dark Energy length scale. While our analysis is done only 
to first order in the perturbation of pu, these features of p}j hold to higher 
orders as well. 

The perturbation, p]j, itself depends on two terms. The first has a power 
law dependence of (r/z/r)^, while the second has a power law dependence of 
(rn/r)^^^- Thus, near the galactic core p]j ~ 1/r^/^, while for r ;» rn, p]i ~ 
1/r^. This behavior at large r is consistent with the pseudo-isothermal density 
profile [P observed in spiral galaxies. Both terms in the perturbation decrease 
rapidly with r, and thus detail structural features of the galaxy disappear 
quickly with distance from its core. 

3.2.2 A Natural Cutoff for p{r) 

With the boundary condition pii{rji) = piii{rii), the only solution to Eq. P^ 
in Region III that is spherically symmetric and finite as r — > c» is 

Pndr)= ^i:(a.)^^(l + 4i+-)^^^^ 

(24) 

In this region, the density decreases to zero exponentially fast, and the profile 
essentially cuts off at \/x/(l + 4:^^"-'^)Xde ', there is therefore a natural cutoff 
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built within the theory. The scale of this cutoff is approximately 0.20 Xde for 
ctA = 3/2, while from Eq. (POJ) . r// w 0.27 Xde', to a good approximation ru w 
\/x/(l + 4i+"^i)A£iB- Regions II and III thus overlap, and the asymptotic 
analysis is self-consistent. As •\/x/(l + 4-'-+"-^)Adb is a scale set by the theory 
as opposed to Eq. (j20p . which is set by an approximation scheme, we take 
r// = ^x/(l + 4i+"-i)A_D_E from now on. 



4 A Cosmological Check 

In this section, we use the recent measurements of the Hubble constant, h, 
by WMAP to set a^- This a a is then used to calculate erg and i?200 using 
the model galaxy derived in the previous section. These calculated values are 
then compared with WMAP measurements for trg, and observations of R2Q0, 
thereby providing a cosmological check of the extended GEOM. 

The density profile for the model galaxy in the previous section was cal- 
culated using rjfj^i, as the background metric in the nonrelativistic and weak 
gravity limits. While this approximation is appropriate at the galactic scale, we 
would not, in general, expect our p(x) to be valid at cosmological scales; at one 
point, we would expect the local curvature of the universe to introduce correc- 
tions that were not taken into account in the original linearization of the met- 
ric, 5/ii/. WMAP and the Supernova Legacy Survey put {2k = — 0.011±o.oi2, 
however, while WMAP and the HST key project set f2K = -0.014±o.oi7 HI- 
Measurements have thus shown that the universe is within experimental error 
of being flat, and indeed, WMAP's determination of ft. = 0.732^0032 was made 
under this assumption. As the density profiles of galaxies do not increase at 
large distances from the galactic core, it is consistent — both internal to the 
perturbative analysis and external with experimental observations — to extend 
our solution in Eq. (P^ to cosmological scales. 

Like cosmological models, the spacetime considered in Sec 2 is divided into 
time-slices. Also like cosmological models, the radial coordinate, r, used in this 
paper is the proper distance between two points on the same time-slice. This 
distance is not what is observed in astronomical observations, however. An 
object separated from an observer by a proper-distance can only be observable 
some time in the future, when the object enters the past lightcone of the 
observer. This is an important distinction that will be used when setting a a. 

Our model from Sec 2 is for a single galaxy, while the universe certainly 
contains more than just one galaxy. That we can apply our single-galaxy model 
to observations made on the universe such as those by WMAP is due to the 
large separations between galaxies, and the universal properties of Pasymp- 
Separations between galaxies are observed to be on the mega-parsec scale, 
while the core size of galaxies, rn, are on the kilo-parsec scale. As we showed 
in Sec 2, the detailed structure of the model galaxy is described by pjj, which 
dies off as l/r^ for large r. Much beyond ~ 6.8r/f, knowledge of the detail 
structure of galaxies is washed out, and only pasymp remains. Thus, the single 
galaxy model in Sec 2 is an adequate model of a galaxy in the universe when it 
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is limited to some region centered on the galaxy, and because of the universal 
properties of Pasymp, the extent of this region can be quite large. Indeed, due to 
the large separation between galaxies, it would be appropriate to extend this 
region to the surface of a sphere of radius, d, where 2d is the proper distance 
from the galaxy to its nearest neighbor galaxy. (If the density is required to 
be everywhere continuous, the shape of this region can be adjusted to ensure 
that this requirement is met.) The profile of our model galaxy is applicable 
within this sphere. Moreover, we are able to describe a distribution of many 
galaxies in the universe by replicating the single-galaxy model; the location 
and parameters of the replicated model need only be changed to reflect the 
properties of each galaxy in the distribution. 



4.1 a A and the Hubble Length, A^ 

The density profile for the model galaxy in Sec 2 decreases exponentially to 
zero beyond a distance ru from the center of the galaxy. As ru does not de- 
pend on the detailed structure of galaxies, and because a precipitous decrease 
in the local density of matter has not been observed in the universe, we argue 
below that we must have d < ru between any two galaxies in our distribution 
of galaxies. This restriction determines r//, which in turn will determine aA- 

Consider a two spiral galaxy model where the two galaxies are separated 
by a distance 2d; galaxy A is located at d, while galaxy B is located at — d. We 
require that d be much larger than the size of either galaxy. We also require 
that while the two galaxies may move relative to each other, this motion takes 
place over a very long time scale, and thus d can approximated as a constant 
independent of time. This is a restrictive requirement; observations show that 
galaxies do orbit one another, and there can be significant relative motion of 
the galaxies. We do not need to construct a detail model of the dynamics of 
galaxies — or even the density profile of them — to set a value for a a-, however. 
All that is required is to determine the maximum separation between galaxies 
within the theory, and then to compare this length with the proper-distance 
between galaxies in cosmological models. This can be done within a model 
where the two galaxies do not move relative to one another. 

As with the single galaxy case, we approximate the motion of stars in both 
galaxies as undergoing circular motion about the center of each galaxy. Since 
the separation between galaxies is much larger than the size of galaxies, the 
motion of stars within one galaxy does not affect the motion of stars within 
the other. Thus, Eq. dTl]) still holds near the core of galaxy A after replacing 
/(x) by /^(x— d). It also holds near the core of galaxy B after replacing /(x) 
by /^(x-|-d). Here, f^ and f^ have the same form as /(x) given by w"^''^'(r), 
but with r replaced by |x ± d|, and with ru and vh replaced by the core size 
and asymptotic rotational velocity of the respective galaxy. 

While we can, at least in principle, follow the same analysis as before to 
determine the density profile of the galaxies — the core of the two galaxies 
are well separated from one another — this detailed analysis is not needed to 
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determine a a- Instead, we note that as in the single galaxy case, spacetime 
can be divided into two domains, 2) and 3- In the first, 23, p » Ade/^tt, 
and asymptotic analysis can still be used to determine the density profile of 
the galaxies as p{x) = pasymp(x) + p}^(x - d) + p}f (x + d). Here, Pasymp 
is a solution of Eq ([T51) . but now in the presence of two galaxies instead of 
one. (While the presence of the galaxies does not explicitly appear in this 
equation, the locations of the galaxies determine the boundary conditions for 
p, and thus establish the underlying geometry as being axially symmetric.) As 
before, pjf{x—d) andp}f (x+d) are the first-order perturbations off pasymp(x) 
near the core of each galaxy; they establish the structure of each. Because of the 
large separation between galaxies, pjf and p}f can be solved independently 
using the equivalent of Eq. (j2ip for the particular Pasymp obtained for the two 
galaxy model, and with the origin shifted to the center of each galaxy. 

In the second, 3, P << ^_D_E/27r, and in this domain, Eq. (IT^ holds. As 
such, the density decreases exponentially to zero, and for most of the region, 
we may take p w 0. Note that this behavior of p is set solely by the form of 
Eq. P^ . and does not depend on either the distribution of the galaxies, or 
the detail structure of them. 

Consider now the behavior of the density profile along the line connecting 
the cores of the two galaxies. If the separation, 2d, between the two galaxies 
is smaller than some distance 2(imax, then for all points along this line p > 
Ade/^tt. The density, p(x), is locally maximum at the core, x = — d, of galaxy 
B, and decreases as the distance, |x-|-d|, from the core increases. Since Eq. (fT^ 
does not depend on the detail structure of either galaxy, the density reaches 
a local minimum, pmin = p(0), at the midpoint between the two galaxies. 
The density then increases as the distance increases until it is again a local 
maximum at the core, x = d, of galaxy A. When d < dmax, Pmin > Ade/^tt, 
and in this case S) is simply connected. 

To estimate the scale of dmax, consider the limit of large d. As the separation 
between galaxies increases, Pmin decreases until at one point pmin << ^db/Stt. 
When this happens, the midpoint between the galaxies is contained in 3, and 
S) is no longer simply connected. Instead, S) — T)a U S)b, where 1)a and S)b 
contains the core of galaxy A and galaxy B, respectively, and 'Da H 1)b = 0. 
The two galaxies are separated by a region in 3 where p decreases exponentially 
fast. As such, for large enough separations the density profile of one galaxy 
will not depend on the presence of the other, and we will be back to the single 
galaxy case. For each galaxy, p is now given by Eq. (j22p . with r measured from 
the center of the galaxy, and the rn and vh set by the properties of this galaxy. 
In this limit, "Da and "Db are spheres centered at the core of their respective 
galaxy, and each has radii r//. Since 2(i,„ax is the maximum separation between 
galaxies for which S) is simply connected, we estimate dmax = fii- 

As astronomical observations have not shown that the density of matter 
decreases precipitously at any point along the separation two galaxies, we 
require that in our model of a distribution of galaxies, the separation between 
galaxies must be less than 2rii. 
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Since 2r// is the maximum possible proper distance between galaxies in the 
model, it must be cosmological in scale. It would be physically appropriate to 
identify this scale with Ae_f/ — the proper distance to the cosmological event 
horizon at present day — for the following reasons. Consider an observer on 
galaxy A, which at the present day, to, is separated from galaxy B by a proper 
distance rjj. If 2rjj < XeH, at a finite time, i, in the future galaxy B will 
pass into the past lightcone of galaxy A, and thus be observed by the observer 
on A. If, on the other hand, 2r// > XeH, then for no time t will galaxy B 
pass within the pass lightcone of galaxy A, and will never be observed by the 
observer on A. More importantly, when 2r// > XeH, no physical property of 
galaxy A can be influenced by galaxy B in the future, and thus the presence 
of galaxy B has no physical significance to the observer on A. Galaxy A and 
galaxy B will be causally disconnected. 

While these are good, physical reasons for identifying 2r// with XeH, it is 
not physically possible to measure XeH directly. All physical measurements are 
made on objects within or on the past lightcone of an observer. Because XeH 
is the proper distance to the event horizon at present day, it cannot, even in 
principle, be measured; it is at best inferred. What is physically measurable, on 
the other hand, is Xh, which is also the radius of the Hubble sphere. Although 
it is shown in [21] that for a ACBM model with (J^Mj-i^yi) = (0.3,0.7) and 
h ~ 0.70 the cosmological event horizon is not equivalent to the Hubble sphere, 
1/H is the time it has taken galaxies to achieve their current separations |14) . 
The Hubble scale is the maximum distance that light has traveled in this time. 
Thus, Xh is the maximum separation between galaxy A and galaxy B from 
which an observer on one could have observed the formation of the other at 
some time in the past; it would be possible for the formation of one galaxy 
can influence the formation of the other. 

There are thus also good physical reasons for identifying 2r// with Xh- As 
the difference between XeH and A^ calculated in [21] is very small at present 
day, and because Xh is a natural length scale in cosmology that is explicitly 
measurable, we identifled 2rH = Xh- Since we can express 



rii 



Sirx 



3/2^(1+41+"-^) 



1/2 

Xh, (25) 



this identification sets the value for aA- Here, i7/i — 0.716±o.o55 from 4 is the 
fractional density of Dark Energy in the universe. The solution to Eq. ([25]) 
gives ua = 1.56±o.iOj which is close to the working value of 3/2 for aA used in 
the previous section. Note that to the lowest significant figure, this aA is the 
ratio of an odd to an even integer as well. 
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4.2 (78, -R20O1 and Observational Data for the Galactic Rotation Curves 

Following [22 and |14| . linear fluctuations in the density within a sphere, 5*8, 
of radius 8h~^ Mpc are characterized by 

j(x) = ^Wzi^, (26) 



ip) 



Ss 



where (• • ■)ss denotes the spatial-average over Ss- Of particular interest is cts, 
the rms fluctuations of 6 within Sg , 

al ^ m^nss (27) 

Since ag is used to study the formation of structure in the universe, usual 

calculations of as involves modeling the evolution of density fluctuations in the 

universe as it expands. These calculations are done in Fourier space, and as is 

usually expressed in terms of a power spectrum, Ak , and a window function, 

W{k) {W{r) in real space), that restricts the fluctuations in the density to 

the region Ss- There are a number of subtleties in this calculation, however, 

and in particular, care must be taken in choosing W{k). Indeed, the simplest 

choice 

1 for r < 8/i~i Mpc 



^''^' ^ for r > 8h-^ Mpc, 

called the top hat window, produces a hard cutoff at r = 8h^^ Mpc, and thus 
can introduce spurious fluctuations at scales k ^ h/8 Mpc^^. 

Our focus is not on the formation of structure during the expansion of the 
universe, however. It is on using as to test the predictions of the extended 
GEOM. As we have an explicit density profile for a model galaxy, we can 
calculate as directly from its definition, Eq. (j27p . using the top- hat window 
function; a power spectrum is not needed in our calculation. 

One estimate of the characteristic length scale for the galaxy-galaxy cor- 
relation function |22j is bh^^; thus there is roughly one galaxy within 5*8. 
Consequently, we calculate as with the approximation that there is a single 
galaxy in 5*8 with its core located at the center of the sphere, and with a spe- 
cific set of values, wj^ and r|^, for the asymptotic velocity and core radius of 
this galaxy. 

As p{r) varies significantly across Ss, we begin by finding the average den- 
sity, 

{p)s = {pH0{rH - r))s + {pasympir)0{r - rH))s + {p}i{r)d{r - rH))s, 

l + UA 



1 + 3q;/i 



Piisympiu8)D{us), (28) 



in Ss- Here, 6{x) is the step function, and ug = 8h ^Mpc/x^^^A_D_E- The 
resultant average density is proportional to Pasymp(i*8), with 

D{us) ^ 1 - yP^ + (^T-:^) C - 1(1 + aA) [^Ccos - i^CsJ ys + 
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-(l + a^)yg/^| 



Here, j/s = rH/{8h ^), and 



xt^cos ~ t^ 



L^si 



rt^si 



J^Cc. 



cos[t/logy8] 
sin [i^ log 2/8] ■ 



(29) 



L^cos — o^ o% 

The parameter 



2aA/(l+aA) 



^^■■- 5<-K^) *""""' '^°' 



C 



-2a./(l+a.) ^ 



H 



(31) 



^(aA)x(ayi) V c 

depends explicitly on the detail structure of the galaxy. As C ^ 4 x 10^'^ for 
w|^ = 200 km/s, D{us) ~ 1, and the asymptotic density dominates {p)s, not 
the detail structure of the galaxy. 

This is not the case for {S{x)'^)s, which involves the integration of p^{r) 
over 5*8. Not surprising, it is now the behavior of the density near the core 
that is relevant. Indeed, we find 



ai = -1 
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1 
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yl + l 
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ys 9/4 + i^2 
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2i^ 
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2j/~ 



2j/ 



-Cc, 



L'si 
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32/8 
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4(l + z/2) 



^cos ^sin I" ■^^^cos^sin 



Ml + iy^)ys 

cos[2i/log2/8] 



cos[z^log2/8] - 
sm[iylogys] j > + 
cos[i/log2/8] + 
sin[i^logy8] j } + 

ZZ/Ocos^sin 



z/(C2„, - C2„) - 2CcosCsi„ sin[2z/logj/8] 



(32) 
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To obtain values for v"^ and rjj , we are guided by the operational definition 
of (78 described in [32]. This definition involves choosing a 5*8 centered at a 
given point in the sky, calculating the average mass in it, shifting this center 
to another point on the sky, repeating the measurement, and continuing until 
all points in the sky is, in principle, covered. The set of all such averages then 
forms an ensemble of such measurements, and the rms fluctuation in the mass 
can then be calculated for this ensemblqj. 

Consistent with this operational definition of as, v'^ and r|j should be 
for a representative spiral galaxy for the universe. Such a galaxy would in 
principle be found through a survey of spiral galaxies, which would result in 
an ensemble of asymptotic rotation velocities and core radii for the galaxies in 
the universe. An average for each parameter can then be taken, and identified 
with vjj- and r|^, which in turn can be used to construct the representative 
galaxy. While such a survey has not yet been done, there exists in the literature 
a large repository of measurements of galactic rotational velocity curves and 
core radii [S] - [13]. Taken as a whole, these 1393 galaxies are reasonably 
random, and is likely representative of the observed universe at large. 

Although there have been a many studies of galactic rotation curves in 
the literature, what is needed here is both the rotational velocity and the 
core radius of galaxies. This requires both a measurement of the velocity as 
a function of the distance from the center of the galaxy, and a fit of the data 
to some model of the velocity curve. To our knowledge, this analysis has been 
done in four places in the literature. (The study [33] gives fits to MOND 
rotation curves, but does not list values for rn-) While each of the data sets 
were obtained with similar physical techniques, there are distinct differences 
in their selection of galaxies, in the exact experimental techniques used, and 
in the models to which the observed rotation curves are fitted. In fact, the 
Hubble constant used by each is often different from one another, and from 
the value of 73.2 km/s/Mpc given by WMAP. The reader is referred to the 
specific papers for details on how these observations were made. 

A number of models are in use in the literature to fit the observed galactic 
rotation curves, and they all require at least two parameters to model observa- 
tions. Four different models of these velocity curves are fitted to the data sets 
we use here, and out of these, two of them can be idealized using w"^°'*'(r). For 
these two models there is a one-to-one correspondence between the parame- 
ters used for the fit with vh and th- Data from these fits are then averaged 
to determine v'^ and r|^. 

A summary of the data sets, and how vh and rn are obtained for each are 
as follows: 



^ Measurements of the average mass in Ss are necessarily done at different times. As 
such, a Markovian assumption must be made that the measurements of the average mass 
over time is equivalent to an ensemble of average masses made at equal times. This same 
assumption must be made when determining the representative galaxy. 
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de Blok et. al. Data Set: De Blok and coworkers made detailed measurements of 
60 LSB galaxies |7^ , and fits of the pseudo- isothermal velocity curve were done 
for 30 of them ^. Later, another set of measurements of 26 LSB galaxies were 
made by de Blok and Bosma [B] , of which 24 are different from the 30 listed in 
[5]. Both the data for the 30 original galaxies, and the 24 subsequent galaxies 
are used here. Although the authors used various models for determining the 
mass-to-light ratio in their measurements, we will use the data that comes 
from the minimum disk model, as this was the one model used for all of the 
galaxies in this set. 

De Blok and coworkers were chiefly concern with modeling the density 
profiles of observed galaxies, and as such, a parameter for the profile, pH, along 
with a parameter for the core radius. Re, were used by them. The asymptotic 
value for their rotational velocity is ^J^ttGphR^- Identifying this expression 
with vh: and v3^c with rn, we are able to extract from the de Blok data 
sets values for vu and th along with their standard error. The authors used 
a value of 75 km/s/Mpc for the Hubble constant. 

CF Data Set: In [T^, Courteau presented observations of the rotational ve- 
locity curves for over 300 northern Sb-Sc UGC galaxies, and determined rn 
for each by fitting the curves to three different models of the velocity, one of 
which, v' {r) = (2wc/7r) arctan (r/r^), is similar to the velocity curve for the 
pseudo-isothermal profile used by de Blok and coworkers. Like the pseudo- 
isothermal curve, v^ (r) can be approximated by the idealized velocity curve 
used here. In the limit r ^ rt, u™ « vc, which sets vc = vh- In the limit 
r ^ Tf, v^'^ sa vc{2r/Trrt), which sets rt — 2r/f/7r. 

A fit of these observations was also made to a velocity curve where the 
steepness of the transition from the hub and the asymptotic velocity curves 
could also be taken into account. The model curve has a form in the r — ^ 
limit that not only disagrees with our idealized profile in one specific case, this 
curve does not fit all the galaxies analyzed by Courteau. We therefore did not 
use data from this fit. 

Values for vc and rt for 351 galaxies was obtained through the VizieR 
service (http://vizier.u-strasbg.fr/viz-bin/VizieR). The great majority of the 
rotation velocity curves were based on single observations of the galaxy; only 
75 of these galaxies were measured multiple times, with the majority of these 
galaxies being observed twice. The data set reposited at VizieR contained 
these multiple measurements, and we have averaged the value of vc and rt 
for the galaxy when multiple measurements were done. The standard error in 
the repeated measurements of a single galaxy can be extremely large; this was 
recognized in [T^]. A value of 70 km/s/Mpc was used for the Hubble constant 
by the author. 

Mathewson et. al. Data Set: In [13', a survey of the velocity curves of 1355 
spiral galaxies in the southern sky was reported. Later, the rotation velocity 
curves for these observations were derived in J18J after folding, deprojecting, 
and smoothing the Mathewson data. Each of these velocity curves are due to 
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a single observation. Courteau performed a fit of Mathewson's observations to 
the v'"''^ curve for 958 of the galaxies in [T2] using a Hubble constant of 70 
km/s/Mpc. The results of Courteau's analysis is reposited in VizieR as well. 
Persic, Salucci, and Stel has proposed the Universal Rotation Curve (URC) 
[15] . [in], which has been used by them to analyze the Mathewson data. In ad- 
dition, this model was the third model used by Courteau in [12] to fit both his 
and Mathewson's data. While the URC asymptotically approaches a constant 
velocity, at small r the URC has a r°-^^ behavior, which is different from the 
pseudo-isothermal curve, the v^ curve, and the idealized velocity curve con- 
sidered here. Although the URC has a different power-law behavior at small 
r, the difference is small enough that it is unknown how erg will change if the 
URC is used in its calculation instead of the v^'^'^^^ curve used here. We leave 
this for future research; for our current purposes, we did not consider fits to 
this velocity curve. 

Rubin et. al. Data Set: In the early 1980s, Rubin and coworkers [5] - [TT] 
presented a detailed study of the rotation curves of 16 Sa, 23 Sb, and 21 
Sc galaxies. This was not a random sampling of such galaxies. Rather, these 
galaxies were deliberately chosen to span a specified range of Sa, Sb, and 
Sc galaxies, and as stated in [10], averaging values of the properties of the 
galaxies in this data set would have little meaning. These measurements can 
contribute to a combined data set of all four measurements, however, and we 
have included them in our analysis. While values for the core radii were not 
given, measurements of the rotational velocity as a function of the distance to 
the center of the galaxy were; we are able to fit this data to the same pseudo- 
isothermal rotation curve used by de Blok, et. al. Results of this fit is given in 
Appendix A. A Hubble constant of 50 km/s/Mpc was used by the authors. 

Wanting to be as unbiased and as inclusive as possible, we have deliberately 
not culled through the data sets to select the cleanest of the rotation curves. 
Nevertheless, we have had to remove the data for 27 galaxies from the data 
sets. A list of these galaxies and the reason why they were removed are given 
in Appendix B, where we have also listed any peculiarities of the four base 
data sets. 

While vh is easily identified for all four data sets, determining th is more 
complicated. For the de Blok et. al. data set, published values of Re was 
first scaled by 75/73.2 to account for differences in the Hubble constant; rn 
is then obtained using rn = V^Rc- The same calculation was made using 
the calculated values of Re from Appendix A, but with 50/73.2 instead 
of 75/73.2 to account for differences in Hubble constants. For the CF and 
Mathewson et. al. data sets, published values of Tt are first scaled by 70/73.2 
to account for differences in Hubble constants, and th is then obtained through 
rn = 7rrf/2. 

The values of v*jj and r*^ are then calculated for three of the four basic 
data sets. As the Rubin et. al. data set was not random, v*jj and r*jj was not 
separately calculated for this data set. These measurements were instead in- 



19 



Data Set 


^H 


r*H 


^200 


0"8 


doBlok et. al. (53) 


119.0±6.8 


3.62±o.33 


210±iio 


0.65±o.ii 


CF (348) 


179.1±2.9 


7.43±o.35 


182±84 


0.92±o.i8 


Mathewson et. al. (935) 


169.5±i.9 


15.19±0.42 


224±50 


0.654±o.093 


Combined (1393) 


172.1±i.6 


11.82±o.30 


206±53 


0.73±o.i2 



Table 1 The v^ (km/s) and r'^ (kpc) is listed for each data set along with the resultant 
^200 (kpc) and erg. The number of data points in each data set is listed in parentheses. 



eluded with data from the other three sets to form a Combined data set, and 
w|j and r|^ were calculated for this set as well. The results of these calcula- 
tions are giving in Table [TJ These values for v'^ and r|^ were then used to 
calculate as using Eq. (|5^ , and the results of this calculation are given in this 
table as well. All the data sets give a value for erg that agree with the WMAP 
value at the 95% confidence level. Moreover, the values for erg for all of these 
data sets — calculated from four sets of observations taken over a thirty- year 
span — agree with one another at the 95% confidence level as well. This is con- 
sistent with our expectation that the sample of 1393 measurements of galactic 
velocity curves used here is reasonably random and unbiased. 

With v^f and r'^, it is possible to calculate R200 by taking in Eq. (^5]) 
ug, — ^ U200 — ^200/x^ -^-DB and {p)s = 200pc where the critical density, 
Pc = l.OOe^QQgg X 10"^^ g/cm^ from [3]. The resultant equation is solved 
numerically for i?20o; the results of this analysis is given in Table [T] as well. 
The range of values for i?200 — with i?200 = 206±53 kpc for the Combined 
data set — is consistent with observations. Notice that like the values for erg, 
the values of i?200 for these data sets also agree with one another at the 95% 
confidence level. 

For both ag and i?200, over 90% of the error is due to the error in a a', 
this is true for all the data sets. The values of erg and i?200 are thus sensitively 
dependent on a a- On the one hand, this sensitivity is consistent with a a being 
a power- law exponent. On the other hand, the excellent agreement between our 
calculated values of as and the WMAP value for erg is all the more compelling 
because of it. 



5 Dynamics of Test Particles in the Model Galaxy 

5.1 The Extended and Gravitational Potentials 

Note that Eq. ([5|) may be written in terms of an effective potential, QJ(x) , as 
x = — VQ3, where 



23(x) = '^(x) + c2 logfH[4 + Sttp/Ade]. 



(33) 



Because of the additional terms from the extension of the GEOM on the 
right hand side of Eq. ([SJ, it is not the gravitational potential, <?(x), that 
determines the dynanncs of massive particles; it is instead 93(x). This is an 
important distinction. Our extension of the GEOM can drastically change the 
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metric, and these changes will have broad implications if it is ^(x) and not 
QJ(x) that determine the dynamics. To see this, we calculate explicitly ^(x) 
in Regions I and II. In Region III, r > r//, and motion in this region is not 
physically relevant. 

Integrating V^^ = AnGp in Region I gives 



${r) = .^(O) + 



.?. . 



JH 



\rH 



(34) 



where ^(0) is an overall integration constant. That this is the usual expression 
for the Newtonian potential can be seen from the relation pH — Svj^/ATrGrjj. 
For Region II, we find 
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(35) 



The \/r term in <?(x) is expected from Newtonian gravity, and is due to the 
boundary conditions for <1> aXr — ru] this is true for the constant terms as well. 
The logarithmic term is due specifically to the source, /(r), as expected. It is a 
long-range potential that extends out to r// , and could potentially explain the 
non-Newtonian interaction observed between galaxies and galactic clusters. 
The p}j{r)r^ term is due to the perturbation of the asymptotic density, and 
contains terms ^ 1/r^/^. It is due to both the boundary terms in p\j and the 
boundary conditions for <!>{r). 

For the c^ term in <?(r), r^pasymp "^ y-SQ/i/Ci+a^)^ which increases as r^-^^ 
for aA = 1.56. This would dominate the dynamics of test particles in the 
galaxy if the extended GEOM depended on <?(x) instead of QJ(x). Indeed, if 
the dynamics were determined by ^(x), the resultant motion for stars in the 
galaxy would not agree with observations whatsoever. Instead, this term and 
the r^p]j term in <P are canceled by the additional density-dependent terms 
in Eq. (1551) . To see this, in Regions I and II p >> Ade/'^t^, and we expand 
Eq. (1331) to give 



QJ(x) = <?(0) 
in Region I, while in Region II, 
^(r) = <l'(0) + ^vj, + vl log 
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DE 



2'""\rH) 2a A \8ttph 



(36) 
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(1 + aA)^ A VH 



c^SiaAW^-""^' - &vl\ - (37) 



2(l + 3aA) V r 



(1 + Sa^) 4ayi I xA|,^ V ^£i£; "^ Ade J J 

where we have used p\i/ Pasymp << 1 and Eq. ([TC]) . The effective potential 
thus increases only logarithmically as r increases. This is expected, and is 
consistent with using the rotational velocity curves in constructing /(r). 

The r^'^^ increase in ^(r), if unchecked, would mean that at some point the 
weak gravity approximation used here would be violated. Indeed, it is doubtful 
if the resultant spacetime will be gravitationally stable, and it would certainly 
be inconsistent with WMAP measurements, which indicate that the universe 
is flat. We found that Region III the density p -^ exponentially fast, however, 
and thus <? — >■ exponentially fast as well when r > rjj . This provides a cutoff 
to <?(r), and prevents the gravitational potential from becoming infinite at 
large r. The growth in ^ ends at r = rjj, and as we find that 'l'{rjj) — ^(0) ~ 
0.02c^, the weak gravity approximation is valid for Regions I and II. 



5.2 Matter Density Measurements Under the Extended GEOM 

That the motion of stars in the galaxy is determined by QJ(x) and not by 
<?(x) has far reaching implications. The local density of matter in the uni- 
verse is determine through observations on the influence that this density has 
on the motion of test particles. These test particles can be divided into two 
classes: massive particles, such as the motion of stars in galaxies, and massless 
particles, such as the motion of photons of various frequencies. 

The extended GEOM affects only the motion of massive test particles. 
As noted above, the motion of such particles is determined by the effective 
potential, 2J. As 2J is in turn determined by p}j, and not the background 
density, Pasymp, observations of the rotational velocity curves of a galaxy us- 
ing the motion of the stars in the galaxy will at best be able to determine 
the perturbation on the background density, p\j, and not Pasymp itself. Since 
Pa-sYxapir) » p]i{i') when r >> th, the majority of the mass in the universe 
cannot be seen with these methods. In particular, the motion of stars in galax- 
ies can only be used to estimate p{r) — /Oasymp; the matter in Pasymp is present, 
but cannot be determined in this way. 

The extended GEOM does not affect the motion of massless particles, 
however. Thus, the trajectory of photons still follow the GEOM, and these 
equations are determined by the local metric, ^^j/. In the nonrelativistic and 
weak gravity limits, 5^^ is determined in turned by the potential, ^(x). Con- 
sequently, measurements of the local density of matter using photons — such as 
through gravitational lensing — results in determining the total density, p, and 
not just the small fraction of it in p]j as is possible using massive particles. 
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6 Concluding Remarks 

Given how sensitive our expression for erg is dependent on a a, that our pre- 
dicted values of erg is within experimental error of its measured value by 
WMAP is a compelling result. This is especially true as the data used in 
calculating as was taken by four different groups over a period of thirty years, 
and for purposes that have no connection whatsoever with our analysis. In the 
absence of a direct experimental measurement of a a-, this agreement between 
the calculated and measured values of as provides a persuasive argument for 
the validity of our extension of the GEOM. 

In [I], we performed a gedanken experiment based on the measurement of 
anomalous accelerations in a dilute gas due to the passage of a sound wave 
with wavenumber k ~ 1 cm~^, and an amplitude that is 10% that of an 
ambient gas that has a density of 10^^* g/cm'^. We were able to establish a 
rough lower bound on aA that for Ade = 7.21 x lO^'^" g/cm'^ is 1.56. This 
is equal to the value of a a found in Sec. 3. While this gedanken experiment 
was a crude estimate of the effects of the extended GEOM and while it is an 
open question as to whether such an experiment would be feasible, the fact 
that the value of aA found is so close to this lower bound raises the possibility 
that direct detection and measurement of aA through terrestrial experiments 
may be possible in the near future. 



A Fitting the Rubin et. al. Data Set 

In nil , measurements of the rotational velocity as a function of radius for 60 Sa, Sb and Sc 
spiral galaxies are given. These measurements allows us to fit the data given to a model of 
the rotational velocity curve. Instead of fitting the data directly to the pseudo-isothermal 
velocity curve, vP~^^°{r), as is done in |5], it is more convenient to fit it to (vP"'^°(r))^. 
Moreover, since what is needed is the asymptotic rotational velocity instead of the density 
parameter for the pseudo-isothermal profile, our fit to (t;P"'^°(r))^ = v^c{r) uses vj{ and 
Rc as the two parameters. Here, 

c{r) = 1 arctan ( — ) ■ (38) 

The variance of the fit is then 



(„p-lSO)2 - J^ , 



JZlivV^y-vUr^)?, (39) 



where {{vn i^n)} is the set TV of rotational velocity verses radius measurements for a 
galaxy. A least squares fit gives 

where (• • •) denotes an average over the data points. While an equation for Rq can also be 
given through minimization of Eq. II39I I. the resultant equation for Rq is given implicitly. 
We find it more useful to substitute Eq. II40II into Eq. II39II . and then use iteration to find 
the Rq that minimizes a^ ■ . Using this value for Rc, vi, is then given by Eq. II40I I. 
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The standard error, a-^jj , in vh , is 
^2 



{«S 



p-lSO,2 



7V4(c(r„)2> 



1 



1 



P-ISO O 



while the standard error, cru in Rq is 



Rh + rl 



ctRc 



Here, 



Rc o", p-iso,„ 



VW^ 



r-l + Rl 



{c(rn?) 



{ri + Rir 

r2 + Rl 



2 / c(r„)r2 
(c(r„)2)\i?2+r2 



(r2+i?2^)2 



-2(c(r„)^> 
((-r^°)'c(r„)> 



ij2 






(r2 + i?2^)2 



1/2 



(41) 



(42) 



(43) 



The results of our fits of the Rubin ct. al. data are tabulated in Table [J] The base data 
from 111! used a Hubble constant of 50 km/s/Mpc, and the results given in the table are for 
this value of the constant. Of the 60 galaxies from llll . NGC 6314 and IC 724 could not be 
fitted to a nonzero Re, while the fit for NGC 2608 resulted in a Re that is less than 0.01 
kpc. 



B Data sets 

For the de Blok et. al. data set, the galaxy F568-3 was analyzed twice; we use the analysis of 
F568-3 given by the authors in |5]. In de Blok and Bosma [6], two of the galaxies, F563-1 and 
U5750, also appeared in [5]; we used the values from [5] for these galaxies in our analysis. 
Finally, the radius for DD0185 from [5] was not determined, and we could not include this 
data point in our analysis. Thus, out of 56 galaxies in this data set, 53 were used. 

For the Rubin et. al. data set, we could not find a nonzero radius for two galaxies, and 
one galaxy had a radius less than 0.01 kpc. As this radius was smaller than the resolution of 
their observations, this data point was not included. A total of 57 galaxies were thus used 
from [TTI . 

For the CF and Mathewson et. al. data sets, the vast majority of the data were based on 
single observations. While we may have greater leeway in removing outliers from this data 
set, even here we were circumspect. First, 75 galaxies in |12| were observed multiple times. 
Of these, the galaxies UGC 7234 and UGC 10096 had listed an asymptotic velocity for one 
of the observations that was opposite from the measured asymptotic velocity for the others. 
We assumed that this was a typographical error, and the sign of the anomalous rotational 
velocity is reversed. Second, five galaxies in the CF and Mathewson et. al. data sets had a 
rn = 0, one galaxy had a radius that was 11-sigma out from the mean, and three galaxies 
had a vu that exceeded 8,000 km/s. These are likely indications that the data was not 
sufficiently accurate to allow for a fit of the velocity curve, and these galaxies were removed 
from the data sets. Finally, given that there are only 1393 galaxies combined in the data 
sets, if a galaxy had a V[j or a rg that was five-sigma or more out from the mean, it was 
removed. In the end, 348 galaxies were used in the CF data set, and 935 galaxies were used 
in the Mathewson et. al. data set. A summary of the data points not used in our analysis is 
given in Table [S] 
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Galaxy 


Rc 


ARc 


VH 


AvH 


Galaxy 


Rc 


ARc 


VH 


AvH 


NGC 1024 


0.27 


0.14 


229.42 


9.77 


NGC 4800 


0.18 


0.06 


171.56 


3.42 


NGC 1357 


0.52 


0.14 


268.19 


16.27 


NGC 7083 


0.89 


0.14 


226.51 


2.27 


NGC 2639 


1.02 


0.34 


337.69 


31.31 


NGC 7171 


2.25 


0.36 


251.35 


6.47 


NGC 2775 


0.40 


0.17 


298.98 


8.90 


NGC 7217 


0.19 


0.10 


275.21 


6.56 


NGC 2844 


0.41 


0.09 


167.50 


18.93 


NGC 7537 


0.80 


0.10 


150.06 


2.35 


NGC 3281 


0.44 


0.05 


211.32 


26.51 


NGC 7606 


1.40 


0.30 


279.29 


4.11 


NGC 3593 


0.16 


0.08 


115.28 


14.01 


UGC 11810 


1.54 


0.38 


193.28 


3.85 


NGC 3898 


0.53 


0.06 


254.76 


28.73 


UGC 12810 


3.22 


0.35 


245.73 


1.47 


NGC 4378 


0.13 


0.06 


307.61 


26.60 


NGC 701 


2.49 


0.58 


188.78 


4.86 


NGC 4419 


0.63 


0.03 


211.55 


2.33 


NGC 753 


0.31 


0.11 


208.50 


3.57 


NGC 4594 


1.65 


0.30 


397.24 


10.15 


NGC 801 


0.79 


0.16 


227.64 


4.06 


NGC 4698 


1.85 


0.47 


284.96 


6.34 


NGC 1035 


1.24 


0.09 


150.62 


1.26 


NGC 4845 


0.11 


0.07 


187.54 


0.07 


NGC 1087 


0.54 


0.10 


131.91 


2.54 


UGC 10205 


2.19 


0.27 


272.34 


4.07 


NGC 1421 


0.54 


0.13 


176.42 


3.94 


NGC 1085 


0.29 


0.05 


307.02 


2.11 


NGC 2715 


1.10 


0.22 


151.47 


2.93 


NGC 1325 


1.80 


0.28 


195.55 


2.67 


NGC 2742 


1.10 


0.16 


181.86 


2.36 


NGC 1353 


0.36 


0.18 


218.48 


8.30 


NGC 2998 


1.08 


0.22 


213.85 


3.22 


NGC 1417 


0.40 


0.05 


278.87 


2.36 


NGC 3495 


3.11 


0.46 


206.75 


3.22 


NGC 1515 


0.06 


0.10 


178.35 


10.03 


NGC 3672 


1.74 


0.24 


208.11 


4.03 


NGC 1620 


1.73 


0.25 


241.62 


3.14 


NGC 4062 


0.79 


0.13 


167.88 


2.65 


NGC 2590 


1.30 


0.54 


255.24 


5.33 


NGC 4321 


0.79 


0.35 


208.24 


5.42 


NGC 2708 


1.91 


0.68 


269.92 


9.45 


NGC 4605 


0.97 


0.32 


112.62 


3.42 


NGC 2815 


1.91 


0.68 


269.92 


9.45 


NGC 4682 


1.17 


0.23 


181.17 


2.97 


NGC 3054 


2.41 


0.56 


259.10 


8.30 


NGC 7541 


0.21 


0.16 


195.04 


5.94 


NGC 3067 


0.76 


0.06 


156.80 


1.22 


NGC 7664 


0.65 


0.14 


196.05 


3.07 


NGC 3145 


0.15 


0.07 


257.00 


4.84 


IC467 


1.64 


0.33 


152.42 


3.26 


NGC 3200 


0.42 


0.09 


266.07 


5.43 


UGC 2885 


0.06 


0.10 


266.22 


5.88 


NGC 3223 


1.35 


0.23 


275.29 


5.51 


UGC 3691 


3.04 


0.33 


229.42 


1.31 


NGC 4448 


0.59 


0.11 


207.02 


1.98 













Table 2 Fitted values of Rq (kpc) and vjj (km/s), and their errors for the Rubin et. al. data 
set. 
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Data Set 


Data Removed 


Reason 


dc Blok ct. al. 


DD0185 


VH 


= oo 


Rubin et. al. 


NGC 6314 


Vh 


= 




IC 724 


TH 


= 




NGC 2608 


vh 


<0.01 


CF 


UGC 6534 


VH 


is 35 a out 




UGC 12543 


VH 


is 11 cr out 


Mathewson et. al. 


ESQ 140-G28 


VH 


> 8, 000 km/s 




ESQ 481-G30 


VH 


> 24, 000 km/s 




ESQ 443-G42 


VH 


> 94, 000 km/s 




ESQ 108-G19 


TH 


= 




ESQ 141-G34 


rn 


= 




ESQ 21-G5 


TH 


is 6 (T out 




ESQ 548-G21 


rn 


is 7 (T out 




NGC 7591 


rn 


= 




ESQ 243-G34 


VH 


is 5 cr out 




ESQ 317-G41 


rn 


= 




ESQ 358-G9 


VH 


is 6 (T out 




ESQ 435-G25 


VH 


is 5 cr out 




ESQ 467-G12 


TH 


= 




ESQ 554-G28 


VH 


is 6 cr out 




ESQ 60-G24 


VH 


is 10 cr out 




ESQ 359-G6 


TH 


is 11 cr out 




ESQ 481-G21 


rn 


is 6 cr out 




UGCA 394 


TH 


is 7 cr out 




ESQ 298-G15 


rn 


is 7 cr out 




ESQ 545-G3 


TH 


is 7 cr out 




ESQ 404-G18 


rn 


is 9 cr out 



Table 3 Listed are the galaxies removed from the data sets used in our analysis along with 
the reason for their removal. 
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